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Abstract 

We pursue our investigation of the non-equihbrium dynamics of the Backgam- 
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temperature. We present an analytical study of the scaling behavior of the local 
correlation and response functions of the density fluctuations of the model, and of 
the associated fluctuation-dissipation ratios, throughout the a-regime of low tem- 
perature and long times. This analysis includes the aging regime, the convergence 
to equilibrium, and the crossover behavior between them. 
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1 Introduction 



There is a long tradition in Statistical Physics to illustrate conceptual problems on ex- 
tremely simplified models . The Ehrenfest dynamical urn model , which was consid- 
ered by Kac as 'probably one of the most instructive models in the whole of Physics' 0, 
is an example of these. It was devised to elucidate the problems raised by the objections 
of Loschmidt and Zermelo against the if-theorem, and in particular to understand the 
fine details of how a system relaxes towards equilibrium. 

The Ehrenfest model is defined as follows. Consider N balls, labeled from 1 to A^, 
distributed between two urns. The numbers of balls in the urns are denoted by A^^i and 
N2 = N — Ni. At each time step a ball is chosen at random (i.e., an integer between 1 
and is chosen at random), and moved from the box where it is to the other box. If 
there was initially an unequal number of balls in each urn, the system will relax to the 
equilibrium state, characterized by the binomial law 



such that each urn contains, on average, an equal number of balls. This model, sub- 
sequently studied by Kohlrausch and Schrodinger was finally fully solved by Kac, 
Siegert and Hess g |, 0. 

The present paper belongs to the same tradition. It is devoted to the analytical study 
of a very simple urn model, which is actually a modern avatar of the Ehrenfest model, the 
so-called Backgammon model This model pertains to the recent stream of research 
where the study of the slow dynamics of complex glassy systems is replaced by that of 
simple models, most of them of mean- field type 0. Its aim was to demonstrate that a 
microscopic model characterized by the absence of energy barriers was nevertheless able 
to exhibit a number of the characteristic features of glassy dynamics. 

The Backgammon model generalizes the Ehrenfest model in two directions. First, 
it is easy to realize that the dynamics of the Ehrenfest model is effectively at infinite 
temperature, because there is no constraint on the move of the drawn ball. Therefore, 
in order to consider the dynamics of the model at finite temperature, one has to define 
an energy and impose e.g. a Metropolis rule for the move of the drawn ball: the move 
is allowed with probability one if the energy decreases or stays unchanged, and with 
probability exp{—i3AE) if it increases by AE. The second generalization consists in 
taking M boxes instead of two, and in going to the thermodynamical limit: A^ ^ 00, 
M 00, with a fixed density p = N/M of balls per box. These two generalizations 
define a class of 'thermodynamical urn models', of which the Backgammon model is a 




(1.1) 



2 



special case. In the Backgammon model the energy is chosen equal to minus the number 
of empty boxes [§]. 

Variants of this model have been considered, where either the choice of energy is 



different, or the a priori statistics of the elementary moves is defined otherwise [10, 11 



12, 13|, or yet other rules are imposed |14]. Mentions of the similarity of the Backgammon 



model to the Ehrenfest model can also be found in refs. ITB 



The expected behavior of the model is easy to grasp. Denoting by A'^^ the occupation 
number of the departure box, determined by the selected ball, and by A^^^ the occupation 
number of the arrival box, drawn at random, four elementary moves are possible: 



Nd>l, Na>0: AE = 
N. = hN„ = 0: AE = 



(diffusion). 



Nd = I, Na > : AE = -1 (condensation), 

Nd>l, Na = 0: AE = +1 (evaporation). (1.2) 

The first two moves, corresponding to diffusion of the balls, are independent of tempera- 
ture. The last two ones are respectively favored at low temperature (condensation), and 
at high temperature (evaporation) 

Therefore, at infinite temperature, since once a ball is drawn it is replaced in a 

box chosen at random, equilibrium is attained rapidly, characterized by a multinomial 

distribution of the balls amongst the boxes, a simple generalization of the result for the 

original Ehrenfest model. When temperature decreases, dynamics slows down because 

the rate of rejection of moves of the drawn ball increases. At equilibrium and at low 

temperature, a few boxes contain a large number of balls. Finally at zero temperature 

the only moves allowed are those for which energy decreases (condensation) or stays 

unchanged (diffusion). The dynamics becomes extremely slow as time passes. Indeed the 

ultimate stage when all boxes but two are empty meets the original Ehrenfest model with 

two boxes, except that now, if one box becomes empty, the dynamics stops. Such an event 

occurs with a probability of order P{N) = [see eq. (|1.1|)], so that the time needed 

to empty one of the two boxes in this ultimate stage is exponentially large in A^, and 

infinite in the thermodynamical limit. It is therefore clear that, at low temperature, after 

a transient regime, two time scales are involved in the process: a fast one, corresponding 

to the equilibration amongst the non-empty boxes, and a slow one, corresponding to 

the occurrence of the rare event of emptying a new box. This picture, though simple, 

accounts for the gist of the model. 

^We use on purpose the same words -diflFusion, condensation and evaporation- as used in the de- 
scription of the elementary moves in the Kawasaki dynamics of an Ising chain |p7|] . There is indeed a 
strong analogy between the latter model and a one-dimensional version of the Backgammon model. 
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A virtue of the Backgammon model is that it allows to test a number of concepts 
which appear in the study of the non-equilibrium dynamics of real systems, though it is 
an extreme simplification of the latter. For instance, slow relaxational dynamics, aging 
properties, exploration of the two-time plane, violation of the fluctuation-dissipation 
theorem [T^, etc. can be investigated in this model, either by numerical means P], in 
a straightforward way, or by analytical methods, which turn out to be much harder |]T^, 
10], ^ 0, 0,|2|. 

The relaxation of energy for the Backgammon model was the subject of previous 



studies |T0|, |20|, |2T|, More recently the relaxation of density has been addressed by 
Franz and Ritort ||2^, in order to investigate possible differences between the behavior of 
this observable and that of the energy. In the present work we give a thorough analytical 
study of this question. The methods used are similar to those introduced in ref. [p3| . 
Our presentation will therefore follow the lines of this reference, yet stay self-contained. 

The outline of this paper is as follows. Sec. 2 defines the quantities studied in this 
work and presents the dynamical equations they obey. In Sec. 3 a formal solution of 
these equations is obtained by the method of characteristics. Sec. 4 makes contact with 
the original Ehrenfest model by considering the simple situation of infinite temperature. 
The two main sections of this work are Sec. 5, devoted to equilibrium properties, and 
especially Sec. 6, which contains a detailed investigation of the non-equilibrium a-regime 
at low temperature. Sec. 7 gives a brief summary and a discussion. 



2 Physical quantities and their dynamical evolution 

Consider a finite system, made of M boxes containing particles. Let Ni{t) be the 
occupation number of box number i at time t, i.e., the number of particles contained in 
that box. We have 

M 

Y^N,{t) = N. (2.1) 

1=1 

The Hamiltonian Ti and action S of the system at inverse temperature (3 read 

s = pn = -pY.^N,fi. (2.2) 

i=l 

Instead of addressing the more ambitious goal of describing the dynamics of the model 
by the configuration C = {Ni, . . . , Nm} involving all the occupation numbers, we restrict 
our study to quantities involving the occupation number of one box only. We consider 
the thermodynamical limit (A^ — > oo, M oo), with a fixed density p = N/M, and 
furthermore restrict the analysis to the homogeneous non-equilibrium initial condition 
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where there is one particle per box. We thus have = M, p = 1, and Ni{t = 0) = 1 for 
i = 1,...,M. 



2.1 Reminder: occupation probabilities and mean energy 

Taking box number 1 as a generic box, we denote by fk{t) the probabihty that it contains 
k particles at time t: 

Ut) = Vmh{Ni{t) = k}. (2.3) 

The dynamical equations obeyed by the occupation probabilities fkit) read [see e.g. 
Appendix A of ref. |2^ for a derivation] 

d/i(t) 2 



dt Ait) 
d/o(t) 



/2(t)+/i(t)/o(t)-2/i(t), 



dt 

with the initial value 
and where 



/i(t)-Mt)/o(t), (2.4) 
/fc(0) = 4,1, (2.5) 



A(t) 



1 + (e-^ - l)/o(t), /i(t) = e'^ + (1 - e-^)/i(t). (2.6) 



Eqs. ( p.4|) describe a birth-death process, i.e., a random walk on the fc-axis with 
position-dependent rates. The rate at which a ball enters a non-empty box (k > 1) is 
equal to 1, while it is equal to fi for an empty box. The first process is temperature- 
independent as it should, since diffusion between non-empty boxes is not constrained 
by temperature. The second process depends on temperature. For instance, at infinite 
temperature, /i is equal to 1, because moving a ball to an empty box is always allowed 
(diffusion and evaporation are both permitted). At zero temperature, /i is equal to /i, 
which expresses the fact that moving a ball to an empty box is allowed if the departure 
box contains one ball only (only diffusion is permitted). The departure rate of a ball from 
a box containing k balls is equal to k/A{t). Hence, there is equilibration of k around 
A(t). At low temperature, A(t) ^ (1 — /o(t))^^, which is equal to the average number of 
balls in a non-empty box. 

To summarize, eqs. (p.4|) describe a random walk along the k-axis, in a confining 
potential centered around A(t), which itself is slowly increasing, whenever the number 



5 



of empty boxes increases |T^. It is indeed intuitively clear that while the equilibration 



inside the potential is fast, the evolution of A(t) is very slow at low temperature. This 
defines two time scales, as mentioned in the Introduction, or two regimes, which are to be 
identified respectively with the (3- and the a-regime, characteristic of glassy relaxation [Q. 

Eqs. ( |2.4|) can be formally rewritten as a Markovian system, of the form 

= Y.^Mt).fmm. (2.7) 

where the matrix A^[/o(t), fiif)] satisfies 

Y.^kAh{t)J\{t)] = Q. (2.8) 

fc>0 

Actually, the reaction terms A(t) and /i(t) involve /o(i) and /i(t), which makes the 
dynamics non-Markovian and the system of equations ( p.4| ) or ( p.7| ) non-linear. 

Eqs. ( |2.4|) or (|2.71 ) ensure the conservation of the moments 



^/,(t) = l, {N,{t)) = Y.kfk{t) = l, (2.9) 

fc>0 fc>l 

respectively corresponding to the total numbers of boxes and of particles. Finally, the 
mean energy per box of the system reads E{t) = — (57Vi(t),o) = —fo{t)- 

2.2 Density correlation function 

The density correlation function c{t, s) measures the correlation between the number 
of particles in box number 1 at times s (waiting time) and t (observation time), with 
0<s<t: 

c{t,s) = {Ni{t)N,{s)) - {Ni{t)){Ni{s)) = {Ni{t)Ni{s)) - 1. (2.10) 
By definition, we have 

{m{t)m{s))= jkPToh{N^{t) = k\N^{s) = j}PToh{m{s) = j} . (2.11) 

j,k>l 



Eq. ( p.lO| ) can thus be rewritten as 



c{t,s)= Y.jkm9i:\t,s)-1, (2.12) 



where 



gi'\t,s)=FToh{mit) = k\N,is)=j}, (2.13) 
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or else as 

c{t,s) = Y.k^k{t,s)-l, (2.14) 

fc>i 

with 

j,{t,s) = Y.Jm9l^\t,s). (2.15) 

i>i 

Since the time variable s plays the role of a parameter in the dynamics, both g^\t,s) 
and 7fc(t, s) obey the dynamical equations 

^il^ = j:^Mt),m]g?{t,s), (2.16) 

d^,it,s) ^ j2Mkim),mMt,s), (2.17) 



dt 



with the initial conditions 

Si) 



gr{s, s) = 7,(s, s) = khis), (2.18) 

implying in particular 

c{s, s) = (iVi(s)2) -l = Y.k' fk{s) - 1, (2.19) 



k>l 



m ae; reement with eq. (|2l0| ). Eqs. (|21^ , ^ ) ensure the conservation of probability, 
i.e., 

E#(t,.) = l, E7.(t,^) = l. (2.20) 

A:>0 fc>0 

2.3 Density response functions 

The local density response function r{t, s) is a measure of the change in the density of 
box number 1 at time t, induced by an infinitesimal change in the conjugate variable, 
i.e., in the present case the local chemical potential acting on the same box, at the earlier 
time s. Assume that box number 1 is subjected to an arbitrary time-dependent chemical 
potential, given by a(t) in reduced units. The perturbed action now reads 

M 

S{t) = -/3^5^,(,),o -a(t)iVi(t). (2.21) 

i=l 

The occupation probabilities of box number 1 now depend on a{t): we denote them 
by fk{t) = Prob{A'^i(t) = k}. In the thermodynamical limit, i.e., to leading order as 
M — > oo, the occupation probabilities of all the other boxes {i = 2, . . . , M) are still 
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given by the fk{t). Indeed the mean-field geometry of the model implies that the mutual 
infiuence of any two distinct boxes scales as 1/M. 

Assuming for simplicity \a{t) \ < /5, the (t) obey the dynamical equations 



dt 

dt 
dt 



{k + l)z._(t)/,VW + - Mt) + ki^-it))fkit) (k > 2 

2u4t)f^it)+fi+it)f^it) - (z/+(t) +/i_(t))/f (t), 

/x_(^)/rw-/^+w/o"W, 



with the initial value 



and with the definitions 



/x+(t) 
/i_(t) 
z/+(t) 
z._(t) 



/i(t)W(-a(t)) + (l-/i(t))e-'^+"W, 
l-/o(t) + /o(t)W(a(t)), 
/i(t) + (l-/i(t))W(-«(t)), 
(l-/o(t))W(a(t)) + /o(t)e-^-"W. 



(2.22) 
(2.23) 



(2.24) 



These equations can be derived along the same lines as in ref. ^3 
In the above expressions, 



W(AS) =min(l,e-^^) 



(2.25) 



is the Metropolis acceptance rate associated with a change of action AS = j3AE. This 
rate is not differentiable at AS* = 0. We have indeed 



dW 



dA^ 



-1, 



dW 



dAS 



0. 



(2.26) 



As a consequence, we are led to define the following two local response functions for the 
density 



Introducing the quantities 



we have 



r^(t,s) 



htit,s) 



5a{s) 



i(s)^0± 



6a{s) 



o(s)^0± 



r^{t,s) = Y.khUt,s). 

k>l 



(2.27) 

(2.28) 
(2.29) 
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The dynamical equations for the s) are obtained by taking the functional deriva- 
tive of eqs. ( p.22|) with respect to a{s) for a{s) — > 0"*^. This yields 

' =Y.^M[foit),m]hfit,s), (2.30) 



df 



with the initial conditions 



2 

ht{s,s) = -— -/2(s) +/x(s)/o(s), 
A(s) 

hoi^^s) = -/^(s)/o(s), 

h^is,s) = -(k + l)e-^fo{s)f,+,is) 

+ (1 - fiis))h-iis) + {ke~^fois) + Ms) - l)Ms) {k > 2), 
/ir(s, s) = -2e-^Us)Us) + fi{s)Us) + (Ms) - 

h^{s,s) = -fi{s)fo{s). (2.31) 
We have 

Y^hUt,s)=0 (2.32) 

fc>0 

at all times t > s. Furthermore, eqs. ( p.31|) imply that both response functions coincide 
at initial times: 

r+(s, s) = r-{s, s) = l + i2e-^ - l)/o(s) - /i(s) + 2(1 - e-^)/o(s)/i(s). (2.33) 



2.4 Fluctuation-dissipation ratios 

The fluctuation-dissipation ratio provides a measure of the violation of the fluctuation- 



dissipation theorem, and thus of the departure of the system from equilibrium g p2| . 
In the present case, we need to define two fluctuation-dissipation ratios 

ds 

These definitions contain no explicit temperature dependence, because a factor of tem- 
perature has been absorbed in the definition of the reduced chemical potential a{t). 

We are thus led to investigate the quantity 

dc{t, s) 



9^ k>i 



E^Cfc(t,s), (2.35) 
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with 



Ck{t,s) 



djk(t, s) 
ds 



(2.36) 



For t > s, the dynamical equations obeyed by these quantities are deduced from eq. ( p.l7|) 
and read 



dt 



Eq. ( |2.31| ) imphes 



d 



k>0 '^^ k>0 



(2.37) 



(2.38) 



Moreover, by integrating eq. ( |2.17| ) for the 7fc(t, s) to first order in 6' = t — s, we obtain 
the initial values 

Ck{s, s) = Y.{k - i)MM[m, fi{t)]m, (2.39) 

£>0 



I.e., 



Ck{s,s) 

Ci(s,s) 
Co(s,s) 



k + 1 



Ais) 
2 



fk+iis) + fk-iis) ik>2), 



A{s) 

-Ms. 



f,{s)+fi{s)fo{s), 



(2.40) 



implying 



dc{t, s) 



ds 



2 - ^ + fo{s) (^(.) - 1) . 



t=s 



Ais) 



(2.41) 



3 Generating functions and integral representations 



It is possible to obtain integral representations of the solutions of eqs. (|2.4| , |2.17| , p. 30 
|2.37D by means of the generating functions 



F{x,t) = Y.hit)x\ G{x,t,s) = Y.jk{t,s)x', 

k>0 k>0 

H^{x,t,s) = 5:/^±(^,s)x^ Z{x,t,s) = J2Ck{t,s)x\ 

k>0 k>0 

These functions obey the partial differential equations 



dF(x,t) 
di 

dG{x, t, s) 
di 



(.-i)(f(m)-J^^-v>w). 

-, N N 1 dG(x, t,s) , 

{x - 1) iG{x, t, s) - ^ ' ^ - Y,{t, s) 



(3.1) 
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dH^{x,t,s) 
di 

dZ{x, t, s) 



(x - 1) H^{x,t,s) - \ ^ ^ J 



A{t) dx 



-Yh±{t,s) 



- 1) f Zix, t, s) - - lc(t, s) 



dt ''^ ^^^^ 

The y-functions which enter these equations read 
Yf{t) = (l-e-^)/o(t) = l-l/A(t), 
Y,{t,s) = (l-e-^)f7o(M) + /o(t)^^-^7o(M)^ 



(3.2) 



dt 



dt 



n±(t,s) = {l-e-^)(^h^{t,s) + fo{t) 



dhUt,s) d/o(t) 



dt 



dt 



h^it,s) 



= (i-e-^)(co(M) + /o(t)^^-^Co(M)). 

The initial conditions for eq. (|3.2| ) are derived from eqs. ( p.5| , p.l8| , |2.31| , p.40D : 

F{x,0) = X, 



(3.3) 



G{x, s, s) 



X ■ 



dF{x, s) 



dx 



H+{x,s,s) = (x-l)f-i-^^^ + e-^/o(.)-/i(.) + 2(l-e-/^)/o(s)/i(s)' 



A(s) dx 



H {x, s, s) = {x — 1) 



X f (1 - hismx, s) + e^^fois)^^^ + (1 - e-^)/o(s)(2/i(s) - 1) 



dx 



1 dF{x, s 



1 + s) - (1 - e-^)/o(s) ((1 - x)Ms) + x) . (3.4) 



^^"^'"'"^ - A(.) dx 
The conservation of moments expressed by eqs. (p^ , p.20| , |2.32| , |2.38| ) imphes 

9F(x,t) 



F(l,t) 



dx 



G(l,t,s) = l, iJ^(l,t,s) = Z(l,t,s) = 0. (3.5) 



Eqs. (p.2|) can be solved by the method of characteristics [see e.g. Appendix B of 
ref. [Ell. We thus obtain 



F{x, t) = (l + (x - l)e-"W) e("-^)^(*'°) + / du)C{x, t, u)Yf{u), 

•J 

G{x,t,s) = e(^-^)^(*'")G(l + (a;-l)e"(^)-"(*\s,s) + f' du IC{x,t,u)Y^{u, s),{3.Q) 



and expressions for H^{x,t,s) and Z{x,t,s) similar to that for G{x,t,s), with the defi- 
nitions 

'■^ du 



r(t) 



A(m) 



(3.7) 
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D{t,u) = / dt;e"(^)-"(*), 

J u 

}C{x,t,u) = (1 -x)e"(")-"W+("-i)^(*'"). 



(3.8) 
(3.9) 



The quantities of interest are given by 



c(t, s) 



dG{x, t, s) 



dx 



-1, r^{t,s) 



dH^{x,t,s) 



x=l 



dx 



x=l 



dc{t,s) dZ{x,t,s) 
ds dx 



x=l 

(3.10) 



Using eqs. (3^) and the corresponding equations for H^{x, t, s) and Z{x, t, s), we obtain 
finally 

c(t,s) = D(t,s)-l + (c(s,s) + l)e"(^)-"W-^*dMe"(")-"(*V-,(M,s), 



r±(t,s) = r±(s,s)e"(^)-"(*) - rdMe"(")'"Wy;±(M,s), 

J s 

OS OS Js 



(3.11) 



Setting a; = in eqs. ( |3.6|) , we get integral representations of the following quantities, 
which will be useful in Sec. 6: 



/o(t) = (l-e--W)e-^(*'°)+ rdw/C(0,t,w)F;(w), 

•J 



lo{t,s) 



-D{t,s) 



G{1- e"(")-"(*\ s, s) + du /C(0, t, u)Y^{u, s) 



with 



]C(0,t,u) = e^i-)-r{t)-D(t,u) ^ 

nil 



du 



(3.12) 
(3.13) 



and similarly 

/lo (t, s 



= e-^^'''^H^ (l - e^(^)-"W, s, s) + du /C(0, t, m)F ± {u, s), 
Co{t,s) = e~^(*'")z(l-e"(")-"(*),s,s) +^*du/C(0,t,u)rc(u,s). 



(3.14) 



4 Infinite-temperature behavior 

At infinite temperature, the Backgammon model, which is a 'thermodynamical' gen- 
eralization of the Ehrenfest model, as discussed in the Introduction, is solvable. The 
solutions given in eqs. (|3.6| , |3.11|) indeed become explicit, because of the simplifications 
e-/3 = A(t) = ^{t) = 1. 



12 



4.1 Reminder: occupation probabilities and mean energy 

Since Yf(t) = 0, we have 

F{x, t) = {l + {x- l)e-*)e("-^)(^-^"*\ (4.1) 
The occupation probabihties thus read 

Mt) = ((1 - e-)^ + ke-^) ^' ^, , (4.2) 

yielding, in particular, the mean energy 

E(t) = -/o(t) = -(l-e-V"-'- (4.3) 



4.2 Correlation and response, fluctuation-dissipation ratios 

Similarly, noticing that Ky(t, s) = Yh±(t, s) = Y^{t, s) = 0, we get 

G{x,t,s) = (l + (a;-l)e^-*) (l + (x-l)e-*(l-e-^))e("-^)(^-^"*\ 
H+{x,t,s) = (x-1) 

X "(e« + {x- l)e-*(e-^ - 1)) e'^-m-e-^yt _ ^(--Dd-e-')- 

H-{x,t,s) = (x - l)e(^-^)(^-^"')-* 



)-s-t+o-''-l 



X 



(l + (a; - l)e-*) e^^"^)^ - (l + (x - l)e^-*) e 



Z{x,t,s) = (x-l)((x-l)e-* + l + e-'^je^"-^^^^-°^^+^-*, (4.4) 
We thus have 

c{t, s) = (l-e-2^)e^-*, r±(t, s) = {l-e-^'+'~'-^)e'-\ = (l+e-'^)e^-*. (4.5) 

Note that both response functions coincide for all times at infinite temperature: r~^{t, s) = 
r~{t,s), whereas the generating functions H^{x,t,s) are different from each other. The 
common value of the fluctuation-dissipation ratios is independent of t, and reads 

1 _ -2s+e~»-l 

-^^(M)= ■ (4.6) 

At equilibrium, i.e., for s ^ 1, the correlation and response functions only depend on the 
time difference 6 = t — s, and they fulfil the fluctuation-dissipation theorem (see Sec. 5): 

ce,{9) = rt^{9) = = e-^ (4.7) 

Equivalently, X='=(t, s) in eq. ( [1.6|) goes to unity as s — > oo. 
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5 Equilibrium properties 



5.1 Reminder: occupation probabilities and mean energy 

At any finite temperature and for long enougli times, tlie system converges towards ther- 
mal equilibrium. The occupation probabilities at equilibrium are given by the stationary 
solution of the dynamical equations 



A —1-1- p~^<=q p/3~-Aoq 
-E., = (/0)e,= " (5.1) 

where Aeq, the equilibrium value of A(t), is related to temperature by 

e'^ = l + (Aeq-l)e^-. (5.2) 

This quantity can be identified with the thermodynamical fugacity of the model [see Ap- 
pendix C of ref. |23]]. Note that the most probable non-zero occupancy is approximately 
equal to Aeq, in agreement with the discussion given below eq. ( ^^ . 

At low temperature, we have 

\n(3 + l In^/J-l 
13 ^ 2/32 



Aeq = (3-\n(3 + —^ — + + 



J, , , 1 , ln/3 , ln-/3-ln/3-l , 

£^eq = -1 + - + H -r^ \ (5.3) 



1^ 

Finally, the generating function of the equilibrium occupation probabilities reads 

, , Aeg - 1 + e^^-i)^-'! , , 

^eq(x) = . (5.4) 



eq 



The convergence towards equilibrium is characterized by a relaxation time teq (de- 
noted as t^q-* in ref. |]2^), which is exponentially divergent at low temperature [see 
eq. m)]. 



5.2 Density correlation and response functions 

At equilibrium, the density correlation function c(t, s) = Ceq{9) and response functions 
r^(t, s) = rfq{0) are stationary: they only depend on the time difference 6 = t ~ s. 

These quantities can be derived from eqs. ( p.ll| ), where all integrals become convo- 
lutions. These equations can therefore be solved in Laplace space. Denoting by Ceq(p) 
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the Laplace transform of Ceq{6), and using similar conventions for the other functions, we 
obtain after some algebra 

rtqiP) = KqiP) = Acq - PCeq(p), (5.5) 

which express that both density response functions coincide at equilibrium, and that they 
obey the identity 

= (5.6) 

In other words, the fluctuation-dissipation ratios X^{6) are identically equal to unity at 
equilibrium, expressing thus the fluctuation-dissipation theorem. 

Furthermore, we obtain the explicit results 

/ A e"'^^'' \ 

Acq - (Acq -1){P+ l/(/0)eq) K{p, A^,) - 
(1 + Meq) (Aeq " (Aeq " 1) (p + l/(/o)eq) KiP, Ae^)) ' 

where 

A J ^ A., I d. .M^. e-.<-) ^ e-. g (5.8) 
is the Laplace transform of the equilibrium kernel 

Kie, Aeq) = exp - A^q (l - e-^/^-)) . (5.9) 



The expression ( |5.7|) is a meromorphic function, whose denominator coincides with 
that found in ref. in the case of the energy correlation and response functions. This 
denominator has an infinite sequence of zeros on the real negative axis, which we denote 
by P = —Pk 1 with A; > 1. The relaxation time of the equilibrium correlation and response 
functions tgq (denoted as tf^ in ref. |2^) is given by the inverse of the smallest one: 

teq = 4) • (5-10) 

Pi 

At low temperature, teq is exponentially divergent (a-relaxation), while all the other 
characteristic times remain microscopic, i.e., of order Agq (/^-relaxation). Indeed, we 
have 



teq - ^^^r^ (A.V-'^^/(A.q) + 1 - A 



A2 

eq 

2e^°<i / 2 \ 2e'3 / 2 In /3 + 1 \ , , 
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with 



n>l '^^^ ^ £>0 

The density correlation and response functions read 



/(A) = / —{e^- _i) = V^^VEt7- (5-12) 



ret(^) = E ^^e-i"^ ce,(^) = E (5.13) 
fc>i fc>i Pk 

where the are the residues of the expression ( p.7[ ) at the poles p = — p^^^ In particular 
the values of these functions at equal times {6 = 0) provide the following sum rules for 
the residues A^: 



<(0) = T.A,= , ce,(0) = E ^ = Ac,. (5.14) 

fc>i fc>i 



The only difference between the expressions ( ^.13 ) of the density correlation and 



response functions at equilibrium and those for the energy |£3[ lies in the values of the 
residues: for the density and for the energy, with 

2 

(5.15) 



A,_ /(Aeq-l)(4')-l/(/o)eq)' 



«fc Aeqpf - 1 

Finally, the analysis presented in this section leads to the following description of the 
behavior of the equilibrium density correlation and response functions at low temperature. 



In the /3-regime (^ ~ 1), expression ( |5.7|) simplifies to 



rtq{p) = Acq - PCeq(p) ^ ^ ^ , , (5.16) 



up to an exponentially small correction. We thus obtain 

g-e/Aeq 

^eq 



<(^) ^ , Ceq(^) ^ Aeq - 1 + C^^/^-. (5.17) 



In the a-regime {6 ~ tgq), the expressions (|5.13| ) are dominated by the first term 



{k = 1). The corresponding residue can be estimated as Ai ^ (Acq — l)/tcq. We 
thus have 

rf^iO) ^ ^^^^e-^/*-, Ceq(e) ^ (Aeq - l)e-^/*-. (5.18) 

In the crossover between both regimes (1^0^ teq), the density correlation and 
response functions exhibit the plateau values 

Ae^- 1 



(4)P1 ^ ^7 ' (Ceq)pl ^ Aeq - 1. (5.19) 



"eq 
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6 Non-equilibrium behavior at low temperature 



In this section, we shall extend to density fluctuations the low-temperature analysis de- 
veloped in refs. I^T], This analytical approach is valid at low temperature, throughout 
the non-equilibrium a-regime (s ^ 1, t — s ^ 1), irrespective of the relative values of s 
or t with respect to the relaxation time teq- 



6.1 Reminder: mean energy 

In the case of the mean energy, our approach consists in simplifying the flrst integral 
equation of eq. ( |3.12|) , by expanding the integrand as a power series in t — u. Keeping 
consistently the flrst two terms of this expansion, we are left with the following differential 
equation 

/o(t) ^ e-^W + Mt)Yfit) + Mt)^^, (6.1) 

with 



roc 

jJt) = / deK{e,A{t)) = K{0,A{t)) = 1 -e 
Jo 

Mt) = - Tie eKie,Mt)) "'-'^^ 
Jo 



dp 



= -A(t)e-^Wj(A(t)). (6.2) 

p=0 



This approach leads to the evolution equation |23 



^^A(A(t),Ae,), (6.3) 



with 

^(A.AJ._L.(i-J^^). (6.4) 

Eq. ( |6.3|) describes the relaxation of energy throughout the a-regime at low temper- 
ature. In particular, 

• At zero temperature, and more generally in the aging regime (1 <^ t <^ teq), we 
have 

dA 1 , , 

hence 

e>o 



dA'/(A')=/ ^(e^^-l-A^) = E^^Vw"T^^^' (6-6) 
Jo Jo ^-^n{n + iy. A^ A^ 



or else 

. / N , , , In In t — 2 , 
A(t) ^Int + lnlntH \ (6.7) 
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In the opposite limit (t ^ tgq); the convergence of A(t) towards its equihbrium 
value Aeq is exponential, 

Aeq-A(t)~e-*/*-, (6.8) 



icq 

with a relaxation time 

A,q_- 1 



^eq ~ . -^(Acq; 



e 



/I \ e^ / 21n/3+ 1 



1 + T^ + --- U + + ■■■ (6-9) 



Acq V K J (^'\ /5 



which is roughly equal to half the relaxation time teq of the correlation and response 
functions, given in eq. ( |5.11| ). 

6.2 Correlation and response, fluctuation-dissipation ratios 



The low-temperature approach consists in simplifying the integral equations ( p.l2| , p.l4| ) 
for 7o(t, s), h^{t, s), and Co{t, s) into the differential equations 



dt 



dt ' 

Co{t,s) ^ j,^t,s)Y^{t,s) + Mt,s)^^^^. (6.10) 
These expressions lead, after some algebra, to the differential evolution equation 

_i_fMM._B(AW.AJ, (6.11, 

where (f){t, s) denotes either 70 (t, s) — fo{t), h^{t, s), or ^o(^) -s), and with 

^<^-^°-)= A(A^eVA) + l-Ar ■ 

Changing time variables from s and t to A(s) and A(t), we obtain the alternative form 

with 
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The normalized solution of the differential equation ( p.l3| ) reads 



/A(s) 

This result can be recast into a multiplicative scaling law 

$(A(5),Aeq) 



A(t) \ 

exp ( - / dA' a(A', Aeq) . (6.15) 



Ht, s) ^ ^),) C A \ , (6.16) 

^ <l>(A(t),Aeq)' ^ ' 



with 



<1>(A, Aeq) = exp dA' a (A', Aeq) 



(6.17) 



Finally, using eqs. ( 3.11| ) and the definitions of the F-functions, we obtain the scaling 
predictions 

'''' (6.18) 



c{t,s) _ r^{t, s) _ (9c(t, s)/ds _ A(t) 'l'(A(s),Ae 
^ rtiis) {dc/ds)^,{s) A{s) <l'(A(t), Aeq) 



throughout the non-equilibrium a-regime, which constitute the main result of this sec- 
tion. The plateau values Cpi(s), 'rpi(s), and {dc/ds)pi{s) in eq. (|6.18|) are the initial 
conditions for the dynamical equations ( |6.11| , |6.13| ). These plateau values generalize the 
results ( |5.19|) to the generic non-equilibrium situation at low temperature. They are at- 



tained when the fast (/3) modes are extinct, i.e., for 1 <^ 6^ = t — s -C s. These quantities, 
which are non-trivial in general, remain to be determined, since their evaluation does not 
pertain to the analysis of the a-regime stricto sensu. 

In the case of the correlation function c{t, s), the plateau value 

Cpi(s) = A{s) - 1 (6.19) 



is simply obtained by replacing Aeq by A(s) in the expression (|5.19|) at equilibrium. More 
generally, we have 

c{t, s) ^ Ceq(t - s; A(s)) ^ A(s) - 1 + e'^^-^^/^^^) (6.20) 



see eq. ^.18|) ] in the non-equilibrium /3-regime. This expression amounts to saying that 



the system is somehow at an instantaneous equilibrium described by A(s). This simple 
description is only justified for moderate times (t — s ^ s), and for quantities, such as 
c(t, s), which remain of order unity for long waiting times and/or at low temperature. 

Figure 1 illustrates the validity of the above description. The zero-temperature corre- 
lation function c{t, s) is plotted against ln(t — s), for different values of the waiting time s, 
indicated on the curves. The full lines show the result of a direct numerical integration 
of eqs. ( |2.4| , |2.17| ). The dashed lines show (to the left, in the /^-regime) eq. ( |6.2(J| ) and 
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(to the right, in the a-regime) the scahng law (|6.18| ), with the plateau value (|6.19| ). The 
quantitative agreement between the exact numerical data and their analytical description 
is increasingly convincing as the waiting time s increases. 

The evaluation of the other plateau values, namely {dc/ds)pi{s) and rpi(s), requires 
more care, since these quantities are already exponentially small at equilibrium. A sketch 
of the method used for their derivation is given in Appendix A. 




Figure 1: Plot of the zero-temperature correlation function c{t,s) against ln(t — s), for 
different values of the waiting time s, indicated on the curves. Full lines: exact numerical 
results. Dashed lines: left (/5-regime): eq. ( |6.2q ); right (a-regime): scaling law ( |6.18| ), 
with plateau value (|6.19|) . 



First, one finds 



A{sye 



2p-A(s) 



A(.) - 1 



.ds)j'^ ^ A(.)2e-AW/(A(.)) + l-A(.) teq[A(.)]' 



(6.21) 



where teq [A(s)] is obtained from eq. ( p.ll| ) by replacing Aeq by A(s). Though the expres- 
sion (|6.21|) is simple, since it can be obtained by replacing Aeq by A(s) in the equilibrium 
result (|5.19|) , its derivation is not. 

The values of t'^i(s), obtained along the same lines, have the more intricate expressions 



2p-A(s) 



A(s)A(A(.), Aeq) (^(A(.) - l)e-^(^)/(A(.)) - 1 + ^ 
A(s)2e-A(«)/(A(s)) + l-A(s) 
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A(s)2e-M«)/(A(s)) + 1 - A(s) 



(6.22) 



The scaling result ( |6.18| ) implies that the fluctuation-dissipation ratios X^{t,s) stay 
constant and equal to their plateau values 



throughout the non-equilibrium a-regime, i.e., for t — s ^ 1. Eqs. (|6.21| , |6.22|) lead to 



.'^.JW 1 - -""-^(^M. A.) - 1 - ^ H- ^) . (6.24) 

These expressions interpolate between the equilibrium values X^(s) = 1 for s teq, and 
the non-trivial values 

+ _ _J_ ( 1 \ e^W 2 2 12 

P^^""^ A(s) ^ V A(s)2y) A(s)/(A(s)) A(s)2 A(s)3 A(s)4 ^ ' ' ' ' 

P'^^^ 1^^ ^ A(I) ~ A(^ ) A{s)I{A{s)) ~ A(^~ A(I)3 ~ A(^^'" 

(6.25) 



at zero temperature, and more generally in the aging regime, i.e., for s <^ teq- 

Figure 2 illustrates the above predictions. The zero-temperature plateau values 
X^(s) are plotted against Ins, for s up to 10^. The full lines show the limit values 
limt_^oo -^^(^, -s), obtained by extrapolating the numerical solutions of eqs. ( |2.4| , p. 17 , 



2.30| , p.37| ). The dashed lines show the analytical predictions (|6.25|) , whose accuracy 



becomes extremely high for the larger values of the waiting time s. 

The above results can be compared with the expression for the zero-temperature 
plateau value of the fluctuation-dissipation ratio associated with energy fluctuations. 



derived in ref. 23 



X--^y(s) ^ Hsms)-1) ^ 2 4 6_ 

A(s)2e-^(^)/(A(s)) + (A(s) - 1)2 A(s)2 A{s f A{s)^ ^' ' 

The expressions ( |6.25| , |6.26D of the three fluctuation-dissipation ratios are very similar. All 
these ratios asymptotically converge to the value 1, characteristic of equilibrium, albeit 
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with large correction terms, which only decay as 1/A(s)^ ~ l/(lns)^, i.e., logarithmically 
slow with the waiting time s. More surprisingly, in the present case of density fluctuations, 
the difference 

A(i) V ~ A(s)/(A(s)) J A(s)2 + A(7")3 ' A{sy 

which is entirely due to the discontinuous slope of the Metropolis acceptance rate [see 
eq. ( p.26|) ], and hence dependent on our choice of a microscopic dynamics for the model, 
exhibits the very same logarithmically slow fall-off in l/(lns)^. 



X+{s) - Xpl(s) ^ ^ ( 1 - 1 ^ T77^ + + T73I + " " " ' (6-27) 




Figure 2: Plot of the zero-temperature plateau values X^i^s) (upper data) and Xpi(s) 
(lower data), against Ins. Full lines: exact numerical results. Dashed lines: analytical 
predictions ( |6.25| ). 

To close up, we give a more detailed account of the dependence on both time scales 
s and t of the scaling law ( |6.18|) : 



At zero temperature, and more generally in the aging regime (s <C teq), we have 

, A2e-^/(A) + l 1117 , , 

^(^-^^ = A(A^e-A/(A) + l-A) "2 + 2A-A^'2A^ + ---' ^''^^^ 

hence 

1/2 1 9 \ 



<|.(A(.),oo) ^ (A(.)e^W) l + -r^ + 



A(s) 4A(s)2 ; 
1 /n , / In In s 4 In In s — 9 \ 
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Therefore, in the aging regime (s, t -C teq), the right-hand-side of the scahng 
law (|6.18|) roughly obeys a square-root behavior, namely 

AW *(A(.),.o) rsV' ^ 1 f 1 1 U . . .V (6.30) 



A(s) $(A(t),oo) \tj V 4 V(lnt)2 (Ins)^ 

In the opposite limit (s ^ tcq), describing the convergence towards equilibrium, 
the scaling function $ blows up exponentially, as 

$(A(s),Aeq)^C(Aeq)e^/*-, (6.31) 

with an s-independent prefactor. In this regime, the scaling law ( |6.18| ) therefore 
exhibits an exponential decay, 

A(t) $(A(.),Ae,) ^ (,.,)/,^^ 

A(5)$(A(t),Ae,) ' ^'-'"^ 

characteristic of equilibrium properties. 



7 Summary and discussion 



The present work, a sequel of refs. fl^, |2T], is devoted to a thorough study of the 



density correlation and response functions and of the associated fluctuation-dissipation 
ratios for the Backgammon model, a mean-field dynamical model introduced in the con- 
text of glassy dynamics We thus provide analytical predictions for the most salient 
findings of Franz and Ritort for the relaxational dynamics of the density fluctuations 
of the model, both at equilibrium and away from equilibrium. 

Let us first summarize our main results. Consider the situation at equilibrium and 
low temperature. The correlations in the Backgammon model possesses a spectrum of 
relaxation times with an exponential separation of time scales, between a /3-regime of 
fast relaxation, involving the microscopic time scale, and a slow a-regime, characterized 
by the relaxation time teq ~ 2e^°''/Aeq ~ 2e^//3^, where Aeq ~ /5 is the fugacity of the 
system. The two-time correlation and response functions are stationary, i.e., they only 
depend on the time difference 9 = t — s. The density correlation function reads [see 
eq. (|;T8|)] 

CeqW- + (Aeq-l)e-^/*-V (7.1) 

/3-regime a-regime 
This expression accounts for the superposition of the fast /^-relaxation (^^ ~ 1 or Aeq), 
and the slow a- relaxation {0 ~ teq)- As a consequence, eq. ( [7. 1|) exhibits a well-defined 
plateau value [see eq. ( |5.19| )] 

(Ceq)pl ^ Aeq - 1 (7.2) 
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throughout the crossover between both regimes, i.e., for 1 ^ ^ t^q. The same conclu- 
sions hold for the response, with 

To^iO) = -^Ceq(^), (7.3) 

which expresses the fluctuation-dissipation theorem. 

One of the most striking outcomes of refs. ||2l| , further extended in the present 
work, is that the exponential separation of fast (/?) and slow (a) modes still holds in a 
generic non-equilibrium situation. This simplifying feature allows an analytical treatment 
of the non-equilibrium a-regime at low temperature. Throughout this regime, we predict 
that the density correlation and response functions obey multiplicative scaling laws of 
the form [see eq. ( |6.18| )] 



c{t,s) r^{t,s) dc{t,s)/ds A(t) $(A(s), Acq) 



cpi(s) r^,{s) {dc/ds)^,{s) A(s)<l>(A(t),Aeq) 



(7.4) 



These scaling laws hold at low temperature (teq ^ 1) and throughout the a-regime 
(s 1, t — s ^ 1), irrespective of the ratios s/teq or (t — s)/teq. When these ratios are 
small, i.e., at zero temperature, and more generally in the aging regime, the rightmost 
side of eq. (|7^ ) assumes an approximate inverse-square- root law (s/t)^/^ [see eq. (|6.3CI|) ]. 
Conversely, when the above ratios are large, i.e., in the regime of convergence towards 
equilibrium at a low but finite temperature, this rightmost side exhibits an exponential 
fall-off, of the form exp(— (t — s)/teq), characteristic of the equilibrium a-relaxation. 

Moreover, in the case of the correlation function, the crossover between the non- 
equilibrium (3- and a-regimes can be described by a simple additive formula, generalizing 
the equilibrium result ( |7.1| ), as 



e(,.)«e~i+(AM-l)M|i^i£M^. (7.5) 

^ A(S) $(A(t),Aeq) 

p-regime ^ ' 

a-regime 

This result, illustrated in Figure 1, amounts to saying that in the non-equilibrium f3- 
regime the system is somehow at an instantaneous equilibrium, described by A(s). The 
validity of such a description is related to the simplicity of the non-equilibrium plateau 
value [see eq. ( |6.19| )] 

Cpi(s) = A(s) - 1. (7.6) 

The same holds for the derivative dc{t, s)/ds. On the contrary, for the response functions 
r^(t, s), there is still a well-defined separation of a fast and a slow component away from 
equilibrium, but this separation is non-trivial, as it involves the plateau values rpj(s), 
which are complicated functions of s, even at zero temperature [see eqs. ( |6.22|) ]. 
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It is worth comparing the results of this paper for the correlation and response of the 
density with those found for the correlation and response of the energy ||2l|, First, as 
underlined in ref. ||2^, the density correlation has an appreciable variation over short time 
scales, i.e., in the /3-regime, while the energy correlation remains essentially constant in 
this regime. Second, in the a-regime, the scaling law (|7.4| ) of the correlation and response 
functions of the density is very similar to that of the energy |]23[. As it turns out, both 
scaling functions only differ by the prefactor A(t)/A(s). Finally, in the present case of 
density fluctuations, because of the choice of a Metropolis dynamics, it is necessary to 
introduce two response functions, r=''(t, s), and accordingly two fluctuation-dissipation 
ratios, X^{t,s). This subtlety was not needed in the case of energy fluctuations. A 
surprising outcome of the present study is that the difference X^i{s) — X~i{s), which 
reflects the choice of the microscopic dynamics of the model, persists in the a-regime, 
where it exhibits a logarithmically slow fall-off [see eq. (|6.27|) ]. 

As can be seen on these results, the Backgammon model is helpful in giving clear- 
cut analytical answers to some recurrent questions in the ongoing investigation of the 
slow dynamics of non-equilibrium systems. This knowledge in turn is of interest to 
see whether one can discriminate between non-equilibrium systems by their relaxational 
behavior, encoded in the correlation and response, or equivalently in the behavior of the 
fluctuation-dissipation ratio. Thus, in the present case, is the Backgammon model more 
akin to a coarsening system, or to a mean-field spin-glass model, or yet to another case? 

In some respects, the quantity A(t), which is omnipresent in our description of the 
dynamics of the model, plays the role of the mean domain size L{t), which is the only 
characteristic length scale in a system undergoing phase ordering |^4|, |2^. Another com- 
mon salient feature of the present model and of coarsening systems is that the scaling 



form ( |7.4|) of the correlation function c{t, s) involves the ratio of a function of A(t) to 
the same function of A(s). On the other hand, as far as the fluctuation-dissipation ratio 
is concerned, there is a maximal discrepancy between the present model at zero tem- 
perature, where both fluctuation-dissipation ratios X'^itjs) go asymptotically to unity 
(with a large logarithmic correction, illustrated in Figure 2), and coarsening systems, for 
which the fluctuation-dissipation ratio is observed to go to zero ||2^, ^7\, (The inter- 
pretation of the latter behavior is that the long-time response comes from the movement 
of the domain walls, and therefore becomes more and more negligible as the domains 



grow ||2^, |28[.) This discrepancy is actually due to the fact that the critical temperature 
of the Backgammon model is = 0. As a consequence, the low-temperature phase, 
where aging persists forever, is reduced to the critical point of the model. More gener- 
ally, as will be presented elsewhere P3, the fluctuation-dissipation ratio exhibits special 
features for systems quenched to their critical point. 
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A Derivation of non-equilibrium plateau values 



In this Appendix, we present the method used in subsection 6.2 for the evaluation of the 
plateau values r^i{s) and {dc/ds)pi{s) in the generic non-equilibrium situation. 

Let us take the example of the latter quantity. First, eq. implies 



ds 



{s) = -Ais)iY^Us). (A.i; 

pi 



Second, (F^)pi(s) can be estimated by performing a Laplace transform on eqs. |3.14| ). 
We thus obtain after some algebra 

yas/pj A(s)''e ^'•*''/(A(s)) + 1 — A(s) 

with 

-A{s)x 



Sc(s)=A(s) dx Z(x,s,s), (A.3) 

JO 1 — X 

where Z{x, s, s) is given in eq. (|3.4| ) in terms of the generating series F{x, s) of the 
occupation probabilities. 

Generalizing the equilibrium results derived in Sec. 5, we anticipate that 5'((s) is 
exponentially small in A(s). As a consequence, it is not legitimate to estimate F{x,s) 
by just replacing Aeq by A(s) in the equilibrium result (^.41). In other words, we first 
have to derive a better evaluation of the occupation probabilities in the a-regime at 
low temperature. This can be done by seeking a solution to the first partial differential 
equation of eq. ( |3.2| ) in the form 

Fix, s) = ^^'^ + e(-^)^(^)A(x, (A.4) 

where the first term represents the occupation probabilities at the instantaneous equilib- 



rium described by A(s) [see eq. (|5.4| )], while the correction term A{x, s) is assumed to be 
small. Eq. (|3.2|) yields a differential equation for A{x,s), which can be integrated as 

^ , , dA(s) /■! , e(i-^)^(") - (1 - y)A(s) - 1 , , 

This expression is proportional to dA(s)/ds, i.e., exponentially small in A(s), as expected. 
The result (|A.5|) has several consequences. 

• Setting X = in eq. we obtain A(0,s) ^ -(dA(s)/ds)/(A(s))/A(s). Re- 

markably enough, we thus recover the evolution equation ( |6.3| ) for A(s) itself. 
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By expanding eq. (|A.4|) to second order in (x — 1), we obtain 



c{s,s) = {N,{sr)-1 = Y.k'U{s)-l = 



k>0 



eqj 



(A.6) 

The second moment of the density fluctuations is thus given by the simple formula 
corresponding to the Poisson law of parameter A (s), namely c(s, s) = (A^i(s)2) — 1 = 
A(s), based on the picture of an instantaneous equilibrium described by A(s), up 
to an exponentially small correction, proportional to dA/ds ~ y4(A(s), Agq), 

Similarly, we can evaluate the initial values of the fluctuation-dissipation ratios 
X='=(s,s). We have indeed r^(s,s) ~ dc{s,s)/ds ~ 1/A(s), and 

= l-e-^/o(^) + ^(/i(^)-c(.,.)) 
fl 1 



2 A{s) 



A(A(.),Aeq), (A.7) 



hence 

X^{s, ^ 1 - (^1 + A(.)A(A(.), Aeq). (A.8) 

The violation of the fluctuation-dissipation theorem at initial times is thus expo- 
nentially small, as it is again proportional to dA(s)/ds. We have in particular 

+ . ^ /I 1 \ 1 A(s)e-'^(^) , 

at zero temperature. 

Finally, coming back to our main goal, i.e., the derivation of the plateau value of 
{dc/ds)pi{s), we insert eqs. (|A.4| , |A.5| ) into the expression ( |3.4| ) for Z{x,s,s), and 
then perform the integral ( |A.3|) . All the integrals brought by the correction term 
A{x,s) can be expressed in terms of /(A(s)) and of elementary functions. We are 
thus left with 

(dc\ A(s)2e-^W 

[ds A(s)2e-^Wj(A(s)) + 1 - A(s) ' ^ ' 

i.e., the result announced in eq. ( |6.21[ ). The results ( |6.22| ) have been derived along 
the very same lines. 
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